Velocity distribution of collapsing starless cores, L694-2 and LI 197 
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ABSTRACT 

In an attempt to understand the dynamics of collapsing starless cores, we have 
conducted a detailed investigation of the velocity fields of two collapsing cores, 
L694-2 and LI 197, with high spatial resolution HON J = 1-0 maps and Monte 
Carlo radiative transfer calculation. It is found that infall motion is most active in 
the middle and outer layers outside the central density-flat region, while both the 
central and outermost parts of the cores are static or exhibit slower motion. Their 
peak velocities are 0.28 km s~^ for L694-2 and 0.20 km s~^ for L1197, which could 
not be found in simple models. These velocity fields are roughly consistent with 
the gravitational collapse models of the isothermal core; however, the velocity 
gradients inside the peak velocity position are steeper than those of the models. 
Our results also show that the density distributions are ~ r~^'^ and ~ ^"^'^ in 
the outer part for L694-2 and L1197, respectively. HCN abundance relative to H2 
is spatially almost constant in L694-2 with a value of 7 x 10~^, while for LI 197, 
it shows a slight inward increase from 1.7 x 10~^ to 3.5 x 10^^. 

Subject headings: ISM: individual (L694-2, L1197) — ISM: kinematics and dy- 
namics — star: formation 



1. Introduction 



Starless cores that are on the verge of or undergoing gravitational collapse provide 
important information on the initial stage of star formation. In order to understand the 
dynamics of such starless cores, it is necessary to obtain information on at least density, 
temperature, and velocity distributions. The density and temperature have bee n determined 
by submillimeter/NIR and NH3 observation wi th a certain degree of accuracy (lEvans et al. 



2OOII : iKandori et al. II2005I : iTafalla et al. 1120041 ) 
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However, the velocity field has not yet been e lucidated. The observations of starless cores 



in the molecular lines of CS, DCO+, and N2H+ flLee et al. Ill999l . l200ll . l2004f ) show that the 



inward motion is extend ed to 0.2 — 0.3 pc w ith a speed of 0.05 — 0.09 km s~^, which is derived 
using a two-layer model (IMyers et al. Ill996l ). The analysis of the CS J = 3-2 and J = 2-1 lines 
with the same model for several cores such as L183, L1521F, L168 9B, L1544, and L694-2, 
suggests that the infall speed appears to increase toward the center ( jLee et al. II2004I ). Since 
the two layer model is rather simple, the derived velocity may be an average value along the 
line of sight and its variation with depth should be considered with caution. Moreover, it 
should be noted that in general, CS is significantly depleted in the central regions of starless 
cores (ITafalla et al. II20021. 120041). On the oth er hands, by applying an improved version of 
the two layer ni o del (jPe Vries fc Myers II2005I ) to the interferometric observations in N2H"*", 
Williams et al. ( 20061 ) found that the infall velocity increases toward the emission peak for 
L694-2 and L1544 . This transition might trace the deep central region of the core because 
the emission of N2H+ primarily comes fro m the central part d ue to its low optical depth and 
centrally condensed matter distribution (ITafalla et al. II2006I ). Therefore, it is highly likely 
that these independent observations do not trace the en tire range of the velocity fields in 



the co llapsi ng cores. Sirn i lar sh ortcomings are evident in ITafalla et al. I (120021 ). iKeto et al. 



( I2004J ). and ISwift et al. I (l2006l ). although they have assumed that the infall speed varies 
monotonically with the radial distance and have performed more sophisticated radiative 
transfer calculations. 

In order to derive a reliable velocity field, by overcoming these problems, we selected 
HCN J = 1-0 hyperfine lines as a probe and observed several candidates of collapsing cores. 
Then, we performed radiative transfer calculations by considering different types of velocity 
fields. An advantage of using the HCN J = 1-0 lines is that three hyperfine lines are well 
separated and they cover a wide range (a factor of 5) of the optical depth. In addition, a 
direct comparison between the integrate d intensity rnap of HCN and that of the N2H"'" reveals 
that HCN is not depleted significantly (jSohn et al. II2004J ). For radiative transfer analyses, 
we used the one- dimensional Monte Carlo co de by considering the line overlap effect due to 



the hyperfine splitting of HCN energy levels (IBernes Ill979l : iGonzalez-Alfonso &: Cernicharo 



I993I : iPark et al. 1119991). Th e refore , the capability of th is cod e is almost t h e sara e as that of 



code used bv ITafalla et al. I (l2002l ). iKeto et al. I (120041 ). and ISwift et al. I (120061 ). The only 



difference is that we explored more complicated velocity fields, by relaxing the condition of 
the monotonic increase or decrease of velocity field, which is an implicit assumption of the 
previous studies. 
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Data 



We analyzed two round-shaped collapsing starless cores, L694-2 and LI 197, located at 
02000 = 19'*41'^04/5, 52000 = 10°57'02."0 and ^2000 = 22'^37'"02/3, ^2000 = 58°57'20."6, 
respectively. The ob servations were pe rformed in September 2003 by using the IRAM 30 m 



at a spacing of 11" (ISohn et al 
for the HCN J 



2004h 



The beam size of the telescope was FWHM = 28" 
1-0 and the velocity resolution of the backend was 0.033 km s^^. Data 



were converted to the Tmb scale using a main beam efficiency of 0.82. We set the frequency 
of the three H CN J = 1-0 hyperfi ne lines according to the Cologne Database for Molecular 
Spectroscopy jMiiller et al. Iboosh as follows: 88.6304157 GHz for F = 0-1, 88.6318473 GHz 
for F = 2-1, and 88.6339360 GHz for F = 1-1. 

In order to perform a comparison with the one- dimensional radiative transfer calcula- 
tions, we assumed that two cores are spherically symmetric. It is justified by the spherically 
symmet ric density distributio n from the observation of dust extinction and emission for 



L694-2 flHarvev et al. ll2003al Jbll The HCN J = 1-0 line profiles are distributed concentri- 
cally and they all exhibit blue asymmetry for L694-2. This is the case for L1197 except for 
red asymmetry in south-west region occupying only < 10% of the observed area, which has a 
negligible contribution to the analysis. The observations will be described in detail in a sep- 
arate paper (Sohn et al. , in preparation). The observed spectra were azimuthally averaged 
at concentric radial annuli with a step of 10". We set the centers of the cores to the peak 
intensity positions of N2H+, which are (0", 0") for L694-2 and (0", 11") for L1197 relative to 
the map centers mentioned above. Since the spectra were at rectangular grid points, they 
were weighted by a factor of exp(— ((i/14")^), where d is the distance in arcseconds between 
the radial distance of the grid point and the specified radius and 14" is one half width of 
telescope beam. The resulting spectra are shown in Fig. [T] and will be used as templates for 
the comparison with synthesized line profiles. 



Model of core 



We synthesized line profiles by using the Monte Carlo radiative transfer model and 
compared them with the observed ones in order to derive the model parameters that best 
explain the o bservations. We adopted distances of 250pc and 400 pc to L694-2 and LI 197, 
respectively (ILee et al. II2OOII ). Velocities relative to the local standard of rest were selected 
as 9.61 km s~^ and —3.16 km s~^ for L694-2 and L1197, respectively. The model cores with a 
radius of 0.15 pc for both L694-2 and LI 197 were composed of 30 concentric shells at regular 
intervals. Additionally, these cores were surrounded by an envelope with an adjustable radius 
(see below). We used 10^ model photons for the Monte Carlo code. 
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The kinetic temperature was assume d to be spatia 



the values derived for L1498 and L1517B (Tafalla et al 



we used the following form adopted by iTafalla et al. I (120041 ). 



ly co nstant at lOK, on the basis of 



2004r) . For the density distribution, 



n(r) 



1 + (r/ro)° 



where is the central density, is the radius of the inner flat region, and a is the asymp- 
totic power index. In dynamic models, decre ases with an increase in the central density 
(IFoster fc Chevaler Ill993l : ICiolek fc Basu Il2000l ). To consider this effect, we defined Vo such 
that it depended on ric, 



2.25Ro 



Ro = a/ ^/AnGfimnnc 



(2) 
(3) 



where Ro is the scale radius of the Bonnor-Ebert sphere, a is the sound speed, /i (= 2.33) 
is the mean molecular weight, and tuh is the hydrogen mass. If a is 2.5 in e quation ([H), it 



repre sents approximately the density distribution of the Bonnor-Ebert sphere (ITafalla et al. 
20041 ). In our model, we tested two cases of sound speed, a = 0.2 and 0.3 km s~^, and finally 
adopted the latter. 

Since the integrated intensity map of the thinnest HCN J = 1-0, F = 0-1 line is similar 
to that of the N?H "|" J = 1-0 lirie, it a ppears that HCN is not significantly depleted for these 
two cores at least (ISohn et al. 1120041 ). In order to take into account any minor depletion or 
enhancement of HCN, we assumed the HCN abundance variation of X{r) = Xo[n{r)/no]^ , 
where /3 is the power index and Hq and Xo are the density and HCN abundance relative 
to H2 near the boundary, respectively. Uo were fixed lastly as 5 x 10^ cm~^ for L694-2 and 
1 X 10^ cm~^ for LI 197, respectively. 

The line profiles are most sensitive to the infall velocity field. We characterized the 
velocity field as a A-shaped field. (The reason why we adopted this shape will be discussed 
in next section.) The innermost region within Ri is static, and the infall velocity increases 
with a linear function of the radial distance toward the peak velocity position (Rm) with a 
maximum magnitude of Vm, and then linearly decreases to the outermost layer. We fixed 
the size of the outermost layer and its velocity to 0.13pc < r < 0.15pc and a constant, 
respectively. In fact, the size of the outermost layer is determined on the basis of the first 
several trials. The velocity of the outermost layer can be directly derived from the line 
profile. The fitting procedure will be described in detail in next section. 

In summary, we have seven free parameters in total : and a for density, Xo and f3 for 
abundance, and Ri, Rm, and Vm for velocity. We made reference to the values of previous 
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studi es for density and abundance ( Harvey et al. 2003a| jbl: iGonzalez-Alfonso fc Cernicharo 
1993h . 



Besides these, a few auxiliary parameters were necessary to fit the fine profiles. The 
e-folding widths of the absorption coefficient profile in the inner static region were specified 
as 0.15 km s~^ and 0.13 km s~^ for L694-2 and L1197, respectively; they are larger than the 
pure thermal width by a factor of ~ 2. In the other layers, these widths were automatically 
adjusted such that the line widths of the model spectra closely resembled the observed ones. 
To explain the hyperfine anomalies of HCN, an envelope with larger microturb ulence width 
and lower density must be introduced (IGonzalez-Alfonso fc Cernicharo Ill993l ). The diffuse 
and turbulent molecular clouds su rrounding the cores traced by CO may be considered as 
the envelope (ITafalla et al. II2006I ). We selected the microturbulence width of 0.5 km s~^ 



and the abundance of 5 x 10~^, which were the same for both cores. The density was set to 
10^ cm~^ for L694-2 and 3 x 10^ cm~^ for L1197. Then, we varied the size of the envelope 
to adjust the relative intensity ratio among the hyperfine lines. It is not necessary for the 
envelope parameters to be unique since various combinations of the density, abundance, and 
size are acceptable if the optical depth is almost the same and the excitation temperature 
is sufficiently low and close to 2.7 K. It should be noted that the envelope is introduced to 
account for only the intensity ratio among the hyperfine lines. 

The line profiles produced by this method were convolved with the telescope beam for 
a direct comparison with the observations. The difference b etween the model line profile 
and the observed one was evaluated by (jZhou et al. Ill993l ). We considered all channels 
of three hyperfine transitions on 5 positions, in the velocity range from —0.7 km s~^ to 
0.7 km s~^. 



4. Results 

In the first phase, we adjusted the density distribution and constant HCN abundance 
without infall motion so that the intensities of F = 0-1 at all locations are similar to those 
observed. Further, we attempted to fit the line profiles with the infall velocity fields of the 
constant or monotonically increasing/decreasing function of the radius and found that they 
cannot reproduce the observed spectra. The infall velocities at every grid points were then 
varied by trial and error. It was immediately evident that the infall motion is dominant in 
the middle and outer layers and other parts of the cores are rather static or exhibit a small 
amount of motion. Therefore, we approximated the velocity field with the A-shaped one. 

The A-shaped velocity field could be qualitatively understood by inspecting the char- 
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acteristic features of the observed line profiles shown in Fig. [T] as follows. On the basis of 
the fact that the position of the self- absorption dip in the thickest F = 2-1 line indicates 
the velocity at the outermost part, we can estimate the infall velocity in that region. These 
infall velocities for L694-2 and L1197 are ~ 0.0 and ~ 0.1 km s~^, respectively. On the other 
hand, the velocity in a region close to the central part can be derived from the line width of 
an optically thin N2H+ line, since it is not depleted significantly and, therefore, i ts emission 



mainly originates from the central part in which the density is the highest. From lLee et al. 



([19991), the FWHM of L694-2 and L1197 are found to be 0.27 km s'^ and 0.28 km s'^, 
respectively. These values can be used to derive the contribution of non-thermal motion 
after removing the thermal width. If this motion is attributed to the systematic inward 
motion, the maximum infall velocity will be 0.12 km s~^ for both L694-2 and L1197. Since 
cores with a pure thermal line width are extremely rare, the systematic motion will be even 
smaller. The line shape also provides information on the velocity field in the middle layer. 
The increasing rate of the brightness from the absorption dip to red part is considerably 
slower than to blue part, particularly for L694-2, thereby resulting in a weaker red peak. 
This suggests that the region of Tj, = 1 remains close to the surface layer in the red part, 
implying that the infall velocity increases steeply toward the deeper layer. Then, it must 
decrease again toward the center to satisfy the constraints of the N2H^ lines. 

In the second phase, in order to derive the velocity field in a more systematic manner, 
we reduced the number of parameters describing the velocity field to three, Ri, Rm, and 
Vjn as mentioned in the previous section. First, we adjusted the parameters of the velocity 
field, while fixing the other density and abundance parameters. Thereafter, the density 
parameters were modified followed by the abundance parameters, while the other parameters 
were unchanged. We repeated this process a few times. 

The spectra of the best fit models are displayed in Fig. [T], where one can note that 
almost all the details are reproduced, particularly for L694-2. The density and velocity 
distribution of the models are shown in Fig. [2l The minimum x^'s are estimated to be 2.7 
and 9.2 for L694-2 and L1197, respectively. The distributions of the around the best fit 
model are displayed in Fig. [3] and Fig. HI As expected, the abundance parameters are not 
completely independent, since both are related to the column density of HCN. There seems 
to be (anti-)correlations among the velocity parameters, too. 

To summarize, the infall velocity begins to increase from near the Ri = 0.035 pc for 
both the cores, and reaches the maximum infall velocity of Vm = 0.28 km s~^ at Rm = 
0.085 pc for L694-2 and Vm = 0.20 km s"^ at R^ = 0.10 pc for L1197 as shown in Fig. H 
The central density and power index of L694-2 are around ric = 1.75 x 10^ cm~^ and a = 2.5, 
while those of L1197 are Uc = 1.75 x 10^ cm~^ and a = 1.5, respectively. Our result of 
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L694-2 is similar to t hat derived from the near-infrared extinction and dust emission maps 
dHarvev et aLlboOSal lbll The abundance is uniform with 7x10 ^ for L694-2 {(3 ~ 0.0), while 
it increases slightly toward the center for L1197 {(3 ~ 0.3) from 1.7 x 10~^ to 3.5 x 10~^. 



5. Discussion 

Since we optimized several parameters cyclically, it is uncertain whether they are a set 
of best fit parameters. This may be true because we did not explore sufficient parameter 
spaces. However, since our optimization is based on the qualitative understanding of line 
profile formation, it is unlikely that a quite different solution will exists. The degree of 
coincidence of the synthesized and the observed line profiles also supports this argument, as 
shown in Fig. [H Although a new solution may be derived, but characteristic features such 
as the A-shaped velocity field will remain unchanged. 

The fact that the velocity distribution is concentrated in the middle and outer layers 
is attributed to the selection of HCN hyperfine lines that can probe a wide range of core 
interiors and to the relaxation of the condition of a constant or monotonic decrease/increase 
in the velocity field. Our model naturally explains the "extended inward motion" observed 
for several starless cores since the infall motion is dominant in the middle and outer layers. 
The size of th e infalling layer is ~ 2/?^ ~ 0.2 pc, which is consistent with the observations 



(IMyers II2005I ). The extended inward motion probed by CS toward L1544 and L1551 appears 



to be ex ceptional in that i t extends over very wide area. This may be due to other types of 



motion (ISwift et al. II2006I ). 



The velocity field that we found does not contradict the observation of I Williams et al. 



(120061 ) ■ in which increasing inward motion toward the center in the central 30" region of L694- 
2 is derived. For this, they used the optically thin N2H"'" line to probe the infall velocity of 

< 0.1 km s~^ and therefore the more dominant motion in the middle and outer layers has 
been missed. In our case, we assumed that the velocity in the central region is zero in order 
to reduce the number of parameters. In fact, we have tested several sets of models whose 
velocity fields have non-zero values in the deeper part, and have found that the velocity of 

< 0.1 km s~^ is acceptable. 

The distribution of the infall motion of our model is similar to that of the ambipolar 



diffus ion model and isothermal core collapse model (ICiolek fc Basu ll2000l : lFoster fc Chevaler 



19931 ) in the sense that the models also have the A-shaped velocity field and their infall 
velocity peaks lie just outside the central density-fiat region. However, the velocity gradient 
of our model inside the velocity peak position is steeper than that predicted by the dynamic 
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models. We can observe such a sudden decrease in the infall motion in the dee per layer in 
the so-called first collapse phase of the model used in iMasunaga et al. I (119981 ) . However, 
this occurs when the central part is very dense and opaque at an evolutionary stage later 
than that of our sample. Some mechanisms resisting against gravit y such as rotation ma y 
work; however, they are easily excluded since rotation is very slow (jWilliams et al. II2006I ). 
A n oscillation of starless cores around equilibrium may be able to explain this, as shown 



m iKeto fc Field I (120051 ). where pressure wave reflected from center combined with inward 
motion due to perturbation from outside results in the steep velocity gradient in the middle. 
However the magnitude of the motion is much less than sound speed. In the case of L694-2 
and LI 197, the inward motion is similar to or greater than the sound speed, suggesting 
that they are undergoing dynamical collapse. The mechanism of a steeper pressure gradient 
should be sought, but it is inappropriate to describe it in greater detail at this stage. 



The motion of the ambipolar diffusion model (jCiolek fc Basu 1 120001 ) is very slow and 
therefore it c annot be used to explai n the velocity field derived in this study. This is because 
the model of ICiolek fc Basu I (120001 ) is developed for L1544 when the infall velocity of this 
core was known to be less than 0.1 km s~^. If the magnetic field is reduced to explain 
the infall motion of the cores with ambipolar diffusion, the model will be supercritical and 
reduced to the usual isothermal collapse model. 



6. Conclusion 



By combining high spatial resolution HCN observation toward L694-2 and LI 197 with 
the Monte Carlo radiative transfer calculation, we have found that the infall velocity is most 
dominant in the middle and outer layers of the cores and the infall velocity can be as large as 
0.3 km s~^. The most important feature of our study is that the velocity field is determined 
with such accuracy and fidelity that it is possible to compare dynamics and observations for 
the first time. An extensive investigation of the velocity distribution in collapsing starless 
cores is in progress (Sohn et al., in preparation). The result will clarify the djTiamic status 
of starless cores in greater detail. 



This study was supported by the Ministry of Science and Technology, Korea, under the 
grant R14-2005-058-01002-0. S.H.L. used the PC cluster facility of Korea Astronomy and 
Space Science Institute. Institute. C.W.L acknowledges support by KOSEF ROl-2003-000- 
10513-0 program. 



-9- 
REFERENCES 

Bernes, C, 1979, A&A, 73, 67 

Ciolek, G.E., Basu, S., 2000, ApJ, 529, 925 

De Vries, C.H., Myers, P.C., 2005, ApJ, 620,800 

Evans, N.J.,II., Rawlings, J.M.C., Shirley, Y.L., Mundy, L.G., 2001, ApJ, 557, 193 

Foster, P.N., Chevalier, R.A., 1993, ApJ, 416, 303 

Gonza lez- Alfonso, E., Cernicharo, J., 1993, A&A, 279, 506 

Harvey, D.W.A., Wilner, D.J., Myers, P.C., Tafalla, M., 2003a, ApJ, 597, 424 

Harvey, D.W.A., Wilner, D.J., Lada, C.J., Myers, P.C., Alves, J.F., 2003b, ApJ, 598, 1112 

Kandori, R., et al. , 2005, AJ, 130, 2166 

Keto, E., Field, G., 2005, ApJ, 635, 1151 

Keto, E., Rybicki, G.B., Bergin, A., Plume, R., 2004, ApJ, 613, 355 
Lee, C.W., Myers, P.C., Tafalla, M., 1999, ApJ, 526, 788 
Lee, C.W., Myers, P.C., Tafalla, M., 2001, ApJS, 136, 703 
Lee, C.W., Myers, P.C., Plume, R., 2004, ApJS, 153, 523 
Masunaga, H., Miyama, S.M., Inutsuka, S., 1998, ApJ, 495, 346 

Miiller, H.S.P., Schloder, F., Stutzki, J., Winnevisser, G., 2005, J.Mol.Struct. 742, 215 
Myers, P.C., 2005, ApJ, 623, 280 

Myers, P.C., Mardones, D., Tafalla, M., WiUiams, J.P., Wilner, D.J., 1996, ApJL, 465, 133 
Park, Y.-S., Kim, J., Minh, Y. C, 1999, ApJ, 520, 223 

Sohn, J., Lee, C.W., Lee, H.M., Park, Y.-S. Myers, P.C., Lee, Y., Tafalla, M., 2004, JKAS, 
37, 261 

Swift, J., Welch, W., Francesco, J.D., Stojimirovic, I., 2006, ApJ, 637, 392 

Tafalla, M., Myers, P.C., Caselh, P., Walmsley, C.M., Comito, C., 2002, ApJ, 569, 815 



Tafalla, M., Myers, P.C., Caselli, P., Walmsley, CM., 2004, A&A, 416, 191 

Tafalla, M., Santiango-Garia, J., Myerys, P.C., Caselli, P., Walmsley, CM., Capsi, A., 2006, 
A&A, 455, 577 

Ward-Thompson, D., Scott, P.P., Hills, R.E., Andre, P., 1994, MNRAS, 268, 276 

Williams, J. P., Lee, CW., Myers, P.C, 2006, ApJ, 636, 952 

Zhou, S., Evans, N.J.IL, Kompc, C, Walmsley, CM., 1993, ApJ, 404, 232 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 11 - 



F=0-1 F=2-l F=l-1 F=0-1 F=2-l F=l-1 




Velocity (km s"') Velocity (km s"') 

Fig. 1. — The comparison of spectra from the best fit model (line) and those from the 
observation (histogram) for L694-2 (left) and L1197 (right) as a function of projected radial 
distance, which is noted in the rightmost panels. Observed spectra are shifted so that their 
LSR velocity is zero. 
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Fig. 2. — The density and infall velocity distributions of the best fit model for L694-2 (left) 
and LI 197 (right). 
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Fig. 3. — Distribution of around the best fit solution with of 2.7 as a function of model 
parameters for L694-2 (see text for the best fit model parameters). Panels a, b, and c are 
about velocity field, and d and e are about density and abundance distribution, respectively. 
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Fig. 4. — Distribution of around the best fit solution with of 9.2 as a function of model 
parameters for LI 197. Others are the same as in Fig. [31 



